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We study theoretically the smectic A to C phase transition in isotropic disordered environments. 
Surprisingly, we find that, as in the clean smectic A to C phase transition, smectic layer fluctuations 
do not affect the nature of the transition, in spite of the fact that they are much stronger in the 
presence of the disorder. As a result, we find that the universality class of the transition is that of 
the "Random field XY model" (RFXY). 
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The effect of quenched disorder on condensed mat- 
ter systems has been widely studied for many years [l|- 
both for practical reasons (since disorder is always 
present in real systems) and fundamental ones. Disorder 
can destroy many types of long ranged order (e.g., ferro- 
magnetic order in systems with quenched random fields 
Q), and it can radically change the critical behavior of 
many phase transitions [3( . 

Such effects have been found in, e.g., 
superconductors @, charge density waves 0, 0], Joseph- 
son junction arrays Q, superfluid helium in aerogelQ, 
and ferromagnetic superconductors [Icj . 

Some of the most novel and dramatic effects of 
quenched disorder are found in liquid crystals confined 
in random porous media (ill. [T3 |. These intriguing sys- 
tems exhibit a variety of exotic "Bragg Glass" phases. 
They also undergo unique types of phase transitions [l3j. 
one of which, the Smectic A to Smectic C (hereafter, AC) 
transition [lj, [l5j|, is the subject of this paper. 

In the high temperature phase (the "A" phase), the 
nematic director h (which points along the axis of align- 
ment of the constituent long molecules that make up the 
smectic material), and the normal TV to the smectic lay- 
ers , are parallel. In the low temperature phase (the U C" 
phase), n and N tilt away from each other. 

The AC transition in clean systems was first consid- 
ered by deGennes[l6|, who showed that, if fluctuations of 
the smectic layers could be neglected, the AC transition 
should be in the universality class of the ferromagnetic 
XY model|l2|. 

The effect of layer fluctuations on this result was con- 
sidered later by Grinstein and Pclcovits [1J] , who showed 
that their effect on the AC transition can, in fact, be ne- 
glected, and that, therefore, the AC transition in clean 
systems is XY-\\\c. 

Unfortunately, for reasons not yet well understood, the 
critical region of the AC transition in clean systems ap- 
pears to be extremely small; result, most exper- 
imental systems exhibit a transition that is accurately 
described by mean-field theory [l^|. As a result, no defini- 
tive experimental test of the above predictions has yet 
been made. 



Recently the nature of the AC transition has been 
studied for a liquid crystal confined in uniaxial fl3j and 
biaxial [l9| disordered environments. In these systems, 
the anisotropy essentially freezes the direction of the 
molecular axes, and the AC transition can be described 
entirely in terms of the configuration of the smectic layers 

EE 01. 

In an isotropic quenched random environment (which 
can be realized most simply by putting the smectic in 
aerogel [20j). which we study in this paper, the problem 
is in many ways more difficult, since now both fluctua- 
tions of the molecular direction and those of the layers 
must be addressed. Indeed, it is not even obvious that 
the two phases between which the transition we wish to 
study occurs even exist in d = 3; the stability of the 
A phase in the presence of even arbitrarily weak disor- 
der remains an open question both theoreticallvfl2j , and 
experimentally [ll|. Presumably, similar issues arise with 
the C phase. 

However, if we assume that both the A and C phases 
are stable, then we are able to completely determine the 
nature of the transition between them. We find that, 
if this stability assumption is correct, the layer fluctu- 
ations do not affect the universality class of this transi- 
tion, which proves to be just that of the random field XY 
modeldJSl. 

This implies a substantial quantitative change in the 
universal critical exponents from their values in the clean 
problem. It is known [2 1| that the coefficients in the 
e = 6 — d expansion for the critical exponents of the 
random field XY model are exactly the same as those 
for the e = 4 — d expansion of the clean (i.e., no random 
field) problem. However, since e = 3 in the physical case 
d = 3 for the random field problem, the e-expansion is 
not quantitatively reliable. It is clear, however, that the 
exponents will be quite different from those for the clean 
XY model, as even the first order in e terms change by 
a factor of 3. 

From a quantitative standpoint, the most useful fea- 
ture of our result is that it connects the exponents of the 
AC transition in an isotropic disordered environment to 
those of a random field XY model, as can be experi- 
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mentally realized in, e.g., anti-ferromagnets with substi- 
tutional disorder [23j. 

The remainder of this paper is devoted to demonstrat- 
ing that the AC transition in the presence of isotropic 
disorder is in the random field XY universality class. 

Our starting model is a modification of the model for 
clean smectics near a Smectic A-Smectic C transition [l4j|. 
the Hamiltonian H = H u + H c + H uc for which consists 
of three parts: 



H u = i J d d r 



K{Vi_uf + B ( d z u - \ 



Vti 



Ai ( V -c) 2 + A 2 fv x cY 



+A 3 



Hue. = - \ d d 



+ Dc 2 + 2vc 4 



gic [ d z u - -\Vu\ 



+ g 2 (V 2 L u)x 
-D(c- V ± it) 2 



where we have defined the direction parallel to the aver- 
aged layer normal in the A phase as the i-axis, and the 
plane perpendicular to z as _L. Here c, which is roughly 
the projection of h onto the layers, is the order param- 
eter for the transition. More precisely, it has only two 
non-zero components (i.e., c z {r) = 0), given by 

ct(f) = [Sij - N^Nj^rijir), i e±, (1) 

where N(r) denotes the unit vector along the layer nor- 
mal, given by N = . Note that all terms in the 

' b J \z-Vu\ 

Hamiltonian are rotation invariant. This is because the 
environment is isotropic and rotating the smectic does 
not cost energy. The pieces H u and H c are, respectively, 
just the elastic energies for smectic layer fluctuations and 
molecular reorientations, while H uc couples u and c. 
The fourth term in H c and the last term in H uc , which 



come from a term D(T) \ N 



% induce the AC tran- 



sition via a sign change in the temperature T-dcpendcnt 
coefficient D(T). For T > T AC , D > 0, and the free 
energy is minimized at c = 0, so the system is in the 
A phase. For T < Tac, D < 0, and the free energy is 
minimized at c ^ 0, so the system is in the C phase. 

Now we include disorder. The aerogel exerts a variety 
of random forces on the molecular axes and the smectic 
layers [T3, [l3| ; the most important of them are[l2|, [l3[ 
the "random tilt fields", which tend to align the local 
molecules and the local normals with the random aerogel 
strands. The contribution of these random effects to the 
free energy can be written as[T3, [H| 



/ 



d d 7 



h(f) ■ V±u + h c (f) ■ c 



(2) 



where the quenched random fields h(f) and h c (r) are 
taken to have Gaussian distributions of zero mean, with 



anisotropic short-ranged correlations: 



h l (r)h J {r') = Mt ] 5 d (r-f') 1 
W^W 7 ) = ^6 d (r- r '), 
ht(7)h j (r>) = A'6±6 d (r-r'), 



(3) 
(4) 
(5) 



The first term in equation ^ has been treated in the 
earlier work[l2j on the smectic A phase in isotropic dis- 
ordered environments, where it leads to strong power-law 
anomalous [Hj. The second term is just the random field 
disorder present in the RFXY model (3. |22||. 

To cope with the quenched disorder we employ the 
replica trick @]. We assume that the free energy of the 
system for a specific realization of the disorder is the 
same as that averaged over many realizations. To cal- 
culate the averaged free energy F = InZ, where Z is 
the partition function, we use the mathematical iden- 
tity In Z = lim 



n^o — n — ■ When calculating Z n , we 
can first compute the average over the random fields 
h(f), whose statistics have been given earlier. Imple- 
menting this procedure gives a replicated Hamiltonian 
H r = HI + H r c + H r uc with the effect of the random 
fields transformed into couplings between n replicated 
fields, with the limit n — > corresponding to the original 
quenched disorder problem: 



1 \ - i 2 
B [ d z u a - - Vti a 



5 / *E 

J a=l 

A f 

+K{V\u a f] -T^f J d d r 22 V ± u a -V ± u p , 
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HZ = 
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A-i(v-c a ) +A 2 (v 



Q = l 



j_ x c Q 



+K 3 



dCn 
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2k R T 
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HZ 



d d r 



E 

.a=l 



gic 2 a ( d z u a - ^|Vw Q | 2 



8c a 



+.g 2 (Viw Q )(V • c a ) + g 3 — ^ ■ d z V±u a 



dz 



+D(c a ■ Vj_w c 



(8) 



If we set u a = 0, the entire Hamiltonian reduces to 
Eq. (J7]), which reduces to the RFXY model if K\ = 
Ki = K 3 . An RG analysis shows that departures from 
this "one constant approximation" (i.e., ^1,2,3 = A) are 
irrelevant 24|; hence, in the absence of the u field, the 
transition is in the RFXY universality class. 

The piece H^ Eq. ^ of H which involves u alone is 
precisely the model for smectics A in isotropic aerogel 
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studied in reference [12|. From the analysis of that ref- 
erence, we know that the critical dimension of Eq. ([5]), 
below which the anharmonic terms in Eq. ([6]) become im- 
portant, is 5. On the other hand, the critical dimension 
of H r c Eq. is well knownjl M, SI to be 6. Because 
of this discrepancy between the two critical dimensions, 
a standard e-expansion study of the entire model Eqs. 
(IHHE])" is impossible. Our solution to this quandry is to 
integrate out only the u a fields pcrturbatively in a mo- 
mentum shell RG approach, which is controlled in an 
e = 5 — rf-cxpansion, to obtain an effective model that 
only involves c a . While unorthodox, this approach is 
very much in the spirit of more conventional RG's: we 
are performing a partial trace over some degrees of free- 
dom to obtain a more tractable Hamiltonian in terms of 
the degrees of freedom remaining after the trace. 

The momentum shell RG procedure consists of trac- 
ing over the short wavelength Fourier modes of u a (r) 
followed by a rescaling of the length. We initially re- 
strict wavevectors to lie in a bounded Brillouin zone 
which we take to be the infinite cylinder <| q± |< A, 
— oo < q z < oo, where A ~ l ja is an ultra-violet cutoff, 
and a is the smectic layer spacing. The displacement field 
u a (f) is separated into high and low wave vector parts 
u a (f) = u^(r) + u>{r), where u>(r) has support in the 
thin wave vector shell Ae~ M <| q± |< A, — oo < q z < oo. 
Here di is a differential parameter di <C 1. We first 
integrate out u^(r). This integration is done pcrturba- 
tively in the anharmonic terms in H Eqns. ([6j)- ([HI) - After 
this pcrturbative step, we anisotropically rescale lengths, 
with r± = f^e e , r z = r' z e uli , so as to restore the UV cut- 
off back to A. This is then followed by rescaling the long 
wave length part of the field with u<(r) = u' a {r')e xl . The 
underlying rotational invariance insures that the graph- 
ical corrections preserve the rotationally invariant oper- 
ator d z u a — i ^V« a j renormalizing it as a whole. It is 

therefore convenient to choose the dimensional rescaling 
that also preserves this operator, which is \ = 2 — uj. 

After this procedure we obtain the following RG 
flow equations to one-loop order, ignoring the term 

D (^c a • ViU a j , since wc arc interested in the critical 
point where D vanishes: 



dB 
~di 
dK 

~dl 

dA 

di 

dgi 

di 
dv 

7c 



d + 3 - 3u - —g J B, 



d-l-h> + —g)K t 



d+1 - uj + —g ) A, 



d+ 1 - u>- — gj £fi, 

d — 1 + ui v. 

128 Bv j 



(9) 
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FIG. 1: Schematic representation of all Feynman graphs that 
renormalize B (top diagram) and gi (bottom diagram). In 
each case the gray box represents all possible ways of con- 
necting the portions of the Feynman graphs shown. These 
are identical for both graphs; as a result, the graphical correc- 
tions to B and gi obey {dB / dl) grap h = (dgi/dl) g raph(gi/B). 
This in turn implies that the anomalous elasticity for gi (see 
text) is identical, up to a multiplicative constant, to that for 
B. 



where g is a dimensionless coupling: 



9 = A 



B 
~K5 



C rf _iA 



d-5 



(14) 



where C'd is the surface area of a d-dimensional sphere 
with radius one divided by (27r) d . 

Note that the graphical corrections inside the paren- 
thesis in Eqs. ([9]) and (fT2]l are the same. This is not 
just an an approximation to one-loop order, but exact to 
arbitrary loop order. This can be easily understood by 
analyzing the structures of the Feynman graphs. In Fig. 
Q] the upper graph summaries all the possible graphical 
corrections to (d z u a ) \ V±u a | 2 ; the lower one does for 
(d z u a )c\. The parts inside the two square boxes are the 
same no matter how complicated they are and how many 
loops they have. 

There are no graphical corrections to (V5_Mo,)(V • c a ), 
which is also exact to arbitrary-loop order. This is be- 
cause both terms have one power of c a while all anhar- 
monic terms have even powers of c Q . Therefore, under 
renormalization both 32,3 flow only as a result of length 
and field rescaling. 

The recursion relations for B, K, and A are identical 
with those found for a smectic A in an isotropic disor- 
dered medium in reference [T2I]. This is also exact to all 
orders, since we have not, in our unusual approach, inte- 
grated out the c fluctuations. This means that all of the 
results obtained in (l2| for the long-wavelength behavior 
of these quantities also hold here. We will also make use 
of this fact later. 

To analyze these flow equations wc introduce an addi- 



tional dimensionless coupling: 173 
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Combining flow 
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FIG. 2: RG flows of the dimensionless couplings g and gs from 
equations (|15|l and (|16[) . All initial models starting to the left 
of the stability limit (73 = 8 flow into the g = g* , g$ — fixed 
point, which therefore controls the AC transition. All models 
starting to the right of the stability limit are unstable. 



Eqs. (j9"l [T3")) with the definitions of g and 53 we find 



di 

dm 
di 



€9 ~32 9 
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93)93, 



(15) 
(16) 



where e = 5 — d. These flow equations have four fixed 
points: g* = or f e,.g| = or 8. The RG flows of 
g and g 3 around these fixed points are illustrated in Fig. 
[2 Note that g% = 8 corresponds to the stability limit of 
the system. Linearizing Eqs. (|15[ I16[) around the only 
stable fixed point g* = 4?e, g^ = 0, we find the graphical 
corrections to v vanish exponentially as £ — ¥ 00. This im- 
plies that integrating out u a only gives a finite correction 
to v, even at arbitrarily long wavelengths. Hence, these 
corrections to v coming from the u a fields do not affect 
the nature of the AC transition. 

During each RG cycle the integration over u> also gen- 
erates terms which do not exist in H r c . The most relevant 
ones are produced in the second cumulant by (9 2 u Q )c 2 
and (V]_u a )(V ■ , 
that the terms generated by 
relevant. 

We'll now show that these terms also do not affect the 
nature of the AC transition. We start with the terms 
generated by (d z u a )c^,: 



E E 9i(k) [k B Tk 2 z G(k)S aP + A(k)k 2 z k 2 ± G 2 (k) 

XC a) j(gi)c a ,j(-gi + k)ci3.i(q2)Cf3a(-q 2 - k) (17) 



Elementary power counting shows 
y±u a ) are less 



' dc a 
< dz 
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disordered environment, we can simply use the results of 
[l2| for the wavevector dependences of these quantities. 
Furthermore, since, as noted earlier, there is an exact re- 
lation between the renormalization of g\ and that of B, 
the wavevector dependence of g\ is identical to that of 
B, up to an overall multiplicative constant. 

Using the just noted connections to the work of fl!| . 
we can simply quote fc-dependences of B, K, A, and g\. 



uvb 



s(fc), 5l (fc) K ^ B ; ? 



k z 
k z 



< ki, 



K(k),A(k) oc 



kl nK ' A , 



< k{, 

>fc < i, 



(18) 



(19) 



9 _ Vb+Vk 
Z 2 



where the anisotropy scaling exponent C 
and t]b,k,A > 0. Another result of [l2| is that the ex- 
ponents rjB,K,A are not fully independent, but connected 
by the exact scaling relation: 
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VB 



(20) 



which is implied b y t he fact that g flows to a nonzero 
stable fixed point [12|. Furthermore, there are certain 
bonds on the values of t\b,k that must be satisfied in 
order for the smectic A phase in an isotropic random en- 
vironment to be stable, which is a prerequisite condition 
for the existence of a sharp smectic A-C transition [l2| 
in such environments. It is only meaningful within these 
bounds to discuss the relevance of the terms in formula 
(fr?]) . These bounds are 



i]k + Vb<2, Vk<1, t]b + 57]k>4:. 



(21) 



The first two bounds come from the requirement of long- 
ranged orientational order and the condition for disloca- 
tions to remain confined, respectively. The third bound 
is obtained by combining 77A > with the exact scaling 
relation (|20)l in d = 3. 

Using expressions (fl^l [!§)) we can write equation (|T7|) 
in a scaling form: 




(.Vb-3 V k)/2 




Xc a ,i(qi)c a ,i(-qi + k)cp,i(q 2 )cp,i(-q2 - k), 



(22) 



where fi,2(x) are scaling functions. Clearly, as k — > the 
replica-diagonal term (i.e., the one which contains 8 a p) 
in (|22[) is irrelevant compared to the quartic (v) term in 
H r c , since its coefficient vanishes like k 7 ^ . 

To decide whether the off-diagonal piece is relevant, we 
treat it as a perturbation and calculate its contributions 
to D: 



where G(k) = l/[B(k)k 2 z + K{k)k\ ]. The fc-dependences 
of B, K, A, and g\ arise due to the the nonzero graph- 
ical corrections in the recursion relations Eqs. (| WT2"1) . 
Because, as mentioned earlier, Eqs. (|9"l fTTj) are identical, 
to all orders, with those for a smectic A in an isotropic 
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It is readily shown that this integral converges for d near 6 
if the exponents t)b,k satisfy the bounds (|2"Tj) . Therefore, 
this off-diagonal piece is also irrelevant. 

Now we discuss the terms generated by (V5_u Q )(V-c Q ), 
which also have a diagonal and an off-diagonal part: 

n 

E E 92 [k B TqiG(q)S aP + &{q)qlG\qj\ x 

a,f) q 

QiQjCa,i(o)cpj(-^). (23) 

Here, unlike g\, has no dependence on q since there 
are no graphical corrections to (V5_m q )(V-c q ). Again we 
can rewrite Eq. (|23j) in a scaling form: 



EE 
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<1± J3 \ ~T ] d af3 + <Zj_ J4 -7 

\1±J \9l 



xqiqjc a ,i{q)cpj{-q), 



(24) 



where fz,i{x) are scaling functions similar to fi,2{x). 
Clearly, both terms are subdominant to the quadratic 
terms in H r c as q — > provided that t\b,k are within the 
stability bounds. 

Therefore, we conclude that integrating out u a only 
gives minor corrections to H T C , which do not affect the 
nature of the transition. Therefore, the universality class 
of the transition is just that of the random field XY 
model, as it would be were the full Hamiltonian just H r c . 

In summary, in this paper we've studied the smectic 
A to C phase transition in isotropic disordered environ- 



ment. Our analysis shows that if the smectic phases are 
stable against fluctuations and unbinding of dislocations, 
the universality class of the transition is that of the "Ran- 
dom Field XY Model" . Surprisingly, in spite of the fact 
that the smectic layer fluctuations are large due to the 
disorder, they have no effect on the nature of the tran- 
sition; that is, if the layers can be frozen by some ex- 
perimental means the universality class of the transition 
still remains the same. During this study we developed 
a "partial renormalization group" strategy which proves 
to be very successful. We expect this strategy to be useful 
in dealing with many problems with anharmonic Hamil- 
tonians which involve multiple fields with different criti- 
cal dimensions. 
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